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We compute the sensitivity to dark matter annihilations for the forthcoming large Cherenkov 
Telescope Array (CTA) in several primary channels and over a range of dark matter masses from 
50 GeV up to 80 TeV. For all channels, we include inverse Compton scattering of by dark 
matter annihilations on the ambient photon background, which yields substantial contributions to 
the overall y-ray flux. We improve the analysis over previous work by: i) implementing a spectral 
and morphological analysis of the y-ray emission; ii) taking into account the most up-to-date cosmic 
ray background obtained from a full CTA Monte Carlo simulation and a description of the diffuse 
astrophysical emission; and Hi) including the systematic uncertainties in the rich observational CTA 
datasets. We find that our spectral and morphological analysis improves the CTA sensitivity by 
roughly a factor 2. For the hadronic channels, CTA will be able to probe thermal dark matter 
candidates over a broad range of masses if the systematic uncertainties in the datasets will be 
controlled better than the percent level. For the leptonic modes, the CTA sensitivity will be well 
below the thermal value of the annihilation cross-section. In this case, even with larger systematics, 
thermal dark matter candidates up to masses of a few TeV will be easily studied. 

PACS numbers: 95.35.-|-d, 95.30.Cq, 98.35.-a, 95.85.Pw 


I. INTRODUCTION 

A substantial body of astrophysical data (such as galaxy rotation curves and weak lensing observations. Large Scale 
Structure surveys and the precise data from CMB measurements mm have revealed the crucial gravitational role of a 
new kind of particle, dubbed as Dark Matter (DM). This provides one of the most compelling hints for beyond Standard 
Model (SM) physics. Nevertheless, in order to identify the microphysics nature of this particle, we need information 
on its mass and interaction characteristics. One common strategy dubbed as indirect detection of DM particles aims 
at detecting the secondary stable SM products of a DM annihilation or decay in our Galaxy and beyond, on top 
of the astrophysical background. For popular DM candidates, belonging to the class of Weakly Interacting Massive 
Particles (WIMPs), these SM products are injected in the realm of investigation of very high energy astrophysics. 
In particular, the current and upcoming ground-based telescopes are very appropriate for probing the properties of 
WIMPs in the mass domain of few 100 GeV up to 100 TeV, by looking at the Very High Energy (VHE) y-rays from 
the annihilation of such particles. 

Over the last several years, in order to identify the nature of DM, the Imaging Atmospheric (or Air) Gerenkov 
Telescopes (lAGTs), such as H.E.S.S. [3], MAGIC |3] and VERITAS [S], have conducted rich observational programs in 
the inner region of the Milky Way halo [SHE], nearby satellite galaxies [SHIH] and galaxy clusters [IMI]. So far, the 
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observations conducted by H.E.S.S. in the inner region of the Galactic halo 150 pc) provide the strongest constraints 
for TeV DM particles [7] . This complements the bounds on the annihilation cross section in the few-tens-of-GeV mass 
range, coming from Fermi-LAT (see e.g. those from the null observations of dwarf galaxies in 7 -rays [22], the Galactic 
Gentre and the diffuse 7 -rays emission measurement at intermediate and high latitudes [2S])) and from 

the effect of generating CMB anisotropies at the recombination epoch and their evolution down to the reionization 
era |30H34]. 

The next-generation of lACTs will be the forthcoming large Gherenkov Telescope Array (CTA) that will surpass 
the overall performances of the present experiments. In particular, several estimates of the CTA sensitivity to DM 
annihilation agree on the fact that there will be a substantial improvement in the flux sensitivity compared to current 
lACTs up to one order of magnitude (see e.g. Refs. [35H39]). 

In this paper we provide a new assessment of the CTA sensitivity to DM annihilations in several primary channels 
(DM DM —)■ e“'"e“, T+r“, 66 , and W'^W~) and over a range of DM masses from 50 GeV up to 80 TeV. We 

improve the analysis over previous works [55H55] in several aspects. More specifically, 

o we include for the first time in the context of the CTA, low energy contributions of the 7 -rays fluxes due to inverse 
Compton Scattering (ICS) on the ambient photon background of the from annihilating DM. We demonstrate 
that this is particularly relevant for the determination of the CTA sensitivity in the leptonic channels, especially 
for the DM DM —>■ e+e” mode. 

o we implement a spectral and morphological analysis of the 7 -rays emissions. We demonstrate that the different 
spectral and spatial behaviours of the DM signal compared to the background can be used to substantially 
improve the CTA sensitivity to DM annihilation. 

o we take into account the most up-to-date Cosmic Ray (CR) background coming from a full CTA Monte Carlo 
simulation and a description of the Galactic Diffuse Emission (GDE) estimated from Fermi-LAT data. In 
particular, in order to maximize the impact of the GDE on the CTA sensitivity, we consider an isotropic GDE 
coming from the inner Galactic halo. With this “extreme” choice, we are clearly overestimating the GDE, 
because it is reasonable to expect that the diffuse 7 -ray contamination decreases as a function of the distance 
from the Galactic Center. As a consequence, our results will be conservative from the point of view of GDE 
uncertainties. 

o we study the impact of the systematics errors in the rich observational CTA datasets. Eor example, possible 
sources of systematic uncertainties arise from observational issues such as different observation zenith angles, 
very high energy emitters and starlight gradients in the field of view, instrumental issues such as broken pixels 
and a non-uniform distribution of PMT quantum efficiency, and performance issues such as 7 -ray and background 
acceptances across the field of view and normalisation between the signal and background regions gsj. 

The rest of the paper is organized as follows: in Sec.|ll|we present the main characteristics of the CTA experiment; in 
Sec. jlllj we provide a short description of the properties of the 7 -rays fluxes coming from DM annihilation, the Region 
of Interest (Rol) relevant for our morphological analysis and the expected number of events in the CTA detector; in 
Sec. jlVj we quantify the CR and the GDE backgrounds; in Sec. [V] we present our analysis methodology. In particular 
the implementation of our the spectral and morphological analysis in the CTA likelihood and the quantification of the 
systematic uncertainties in the CTA sensitivity; in Sec. jVIj we present our results and in Sec. m we conclude. 


II. THE CHERENKOV TELESCOPE ARRAY 

CTA will be the next-generation array of lACTs in VHE 7 -ray astronomy. It is envisaged as a two-site array to 
allow full coverage of the sky. One in the northern hemisphere will have the aim of studying extragalactic sources, 
while the other in the Southern hemisphere will emphasize observations of the Galactic Center (GC) region since this 
can be observed close to the zenith during the austral winter. 

CTA will consist of several tens of two or three different types of telescope displaced over a km square area, with 
sizes of about 6 , 12 and 23 m in diameter, respectively. The sensitivity is expected to be a factor 10 better than 
currently operating lACTs: the field of view (EoV) of the small size telescope will be around 9°, the angular resolution 
roughly an arcminute, the energy resolution from about 20% at 100 GeV to better than 5% at 10 TeV, and a lower 
energy threshold of several tens of GeV mi- Although the final design of the array is not settled yet, a detailed 
Monte Carlo study has been performed on various candidate array configurations to estimate the array performances 
in terms of background rejection, point-source sensitivity, angular and energy resolutions |42j . 
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In this study, we will consider the instrument response functions^ obtained for the proposed array 1. This benchmark 
array is a balanced choice to allow for good sensitivity both in low and high energy regime. For this array, we make 
use of the information of the effective area for photons, residual background rate, angular and energy resolutions 
provided in Ref. [32] • 

Observations by lACTs in high energy y-ray astronomy traditionally employ two regions on the sky expected to 
have roughly the same astrophysical emission, but significantly different amounts of DM annihilation. The signal 
region, where the larger annihilation signal is expected to be, is usually dubbed as the ON region. The other, called 
the OFF region, is taken to be larger than the ON region and is used for the background determination. The two 
regions are chosen in nearby regions of the sky and the statistical analysis is performed by using a test statistic 
defined as the difference in counts from the two regions. Observations usually cover part of the ON and OFF regions 
simultaneously in order to avoid additional systematics that may arise from different atmospheric, instrumental or 
observing conditions. 

In this study, we will optimize the above mentioned ON-OFF method by carrying out a full likelihood analysis which 
uses the expected spectral and spatial distribution of the DM signals. As will become clear later on, this improved 
analysis takes full advantage of the spatial morphology of the signal with respect to the background by using the 
available information over several degrees. Nevertheless, due to the limited CTA FoV (around 7.5° for middle-size 
telescopes^), several observational pointings are needed to accurately map the inner 10° of the Galactic halo. In 
particular, due to the radial dependence of the acceptance in the FoV, one needs to develop a proper observation scan 
strategy to reduce the systematics that may arise at the edges of the FoV of the observation in order to uniformly 
map the Rol. 

In what follows, we will consider an optimized observational strategy which provides, in addition to a uniform 
exposure of 500 hours over all the Rols relevant for our analysis, a substantial exposure in Rols, used as OFF regions 
for background measurements, beyond the central 5 degrees. The choice of the pointing strategy to uniformly map 
the considered Rols crucially depends on the total observation time, number of pointings and grid spacing. 

A general and realistic optimization of the observation strategy to search for DM in the GC region is beyond the 
scope of the present paper. For example, the OFF regions for core profiles (e.g. isothermal or burkert), need to be 
taken far away from the ON region in order to enable a significant gradient between the two regions. Given the fact 
that the CTA FoV is limited, dedicated OFF observations are required to estimate the residual y-ray background in 
the DM core. Such observational strategy and analysis methodology are not described here. Therefore it is worth 
stressing that our study, like those presented in Refs. |35H39|, is only suitable for cuspy DM prohles. 


III. DARK MATTER SIGNALS 
A. Properties of the DM fluxes 


Annihilating DM particles induce high energy y-rays fluxes both by direct emission (prompt) and by ICSs of 
produced by DM annihilation on the ambient photon background (secondary). 

• prompt emission: The differential y-ray flux, produced by the prompt annihilation of self-conjugate DM particles 
of mass rriDM, coming from a given angular direction dD, is written as 




dDdi?..v 


1^0 pI 

2 Ltt 


JV»T. 
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dN-f 


7 


■'("> = / 0 

Jl.o.s. ^0 


ds p'^(r(s, 9)) 


pI 


( 1 ) 


where (av) / and dN^/dEj are respectively the self-annihilation cross-section and the energy spectrum of photons 
per one annihilation in the channel with final state /. Here the coordinate r is written r(s,9) = (rQ + s'^ — 
2rQScos0)^/^, where s is the parametrization for the distance along the line-of-sight (l.o.s.), 9 is the aperture 
between the direction of observation and the Galactic plane and tq = 8.33 kpc is the Sun’s location with respect 
to the GC. As usual, the aperture 9 in polar coordinates can be expressed in terms of galactic latitude b and 
longitude £ via the relation 9(b,£) = arccos (cos 6 cos ^). The function J(9), commonly referred as to the J-factor, 


^ Notice that in Ref. m the instrument response functions are taken from unpublished materials. Therefore a direct comparison with 
the results presented here and in previous CTA sensitivity studies is in principle not possible. 

^ The small and large-size telescopes are expected to have a FoV of 9° and 4.5°, respectively. 
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i-esime Rol 

Angular size (ADi) in [sr] 

J-factor plfro) in [GeV^/cm®] 



Einasto used here 

Einasto as in [3^ 

NEW as in gS] 

first Rol (9i = 1°) 

5.91 X 10“'‘ 

1.42 X 10®^ 

4.61 X 10®^ 

1.09 X 10^^ 

second Rol (0; = 2°) 

2.51 X 10“® 

3.17 X 10®^ 

9.56 X 10®^ 

2.02 X 10®^ 

third Rol (Oi = 3°) 

4.42 X 10“® 

3.37 X 10^^ 

9.59 X 10®^ 

2.03 X 10^^ 

fourth Rol {Oi = 4°) 

6.33 X 10“® 

3.30 X 10^^ 

9.02 X 10®^ 

1.98 X 10^^ 

fifth Rol {6i = 5°) 

8.23 X 10“® 

3.15 X 10®^ 

8.35 X 10®^ 

1.92 X 10®^ 


TABLE I: Angular size of the i-esime Rol and corresponding value of the J-factor in units of GeV^/cm® for the Einsato profile 
considered here. For a sake of comparison, averaged J-factors are also given for an alternative normalisation of the Einasto 
profile [39] and a NEW profile [45]. 


integrates the square of the DM density p along the line of sight. For the computation of the J-factor, we will 
always assume an Einasto profile 


p{r) = ps exp 




( 2 ) 


whose parameters {ps = 0.033 GeV/cm^,rs = 28.44 kpc, as = 0.17) are taken by demanding that: i) the 
DM density at the Sun location is pq = 0.3 GeV/cm^; ii) the total DM mass contained in 60 kpc is Mgo = 
4.7 X IO^^Mq [44]. In order to derive the prompt y-rays fluxes, we will compute J{9) and dN^/dEj, by using 
the tools in Ref. |i5] . 

It is worth stressing here that either other parametrizations or parameters (e.g. the authors of Ref. |39| quote 
Pq = 0.4 GeV/cm^ at Tq = 8.5 kpc) of the Einasto profile can yield a substantially different value of the 
J-factor. Therefore, when comparing different results, we have always to be aware about this fine distinction. 
In particular, our present choice of parameters leads to the smallest value of the J-factor in the inner Galactic 
halo compared to previous analyses |35H39]. As a consequence, our results, from the point of view of DM profile 
uncertainties, are the most conservative. 


• secondary emission: The differential y-ray flux produced by the IG radiative processes within an angular region 
dD of the sky can be obtained by convolving the number density in the emitting medium with the differential 
power that it radiates. We write 


dDdi?..y 


1 



dEe'Pic(E^,Ee,r) 


dne± 


{Ee,r) , 


(3) 


where the factor 2 takes into account that an equal population of electrons and positrons is produced by DM 
annihilations. In Eq. dne±/dEf. is the number density of electrons (or positrons) after energy losses and 
diffusion which obeys the diffusion-loss equation given for instance in Refs. [45||46], while Vic = 
total differential power radiated into photons by the IGS mechanism®. We refer the reader to Refs. |151I3H] and 
references therein where the analytic formula of Vic, valid in the full Klein-Nishina case, and a semi-analytic 
solution of drZgi /dEg are provided. 

Eollowing Ref. |3S|, Eq. ([^ can be rewritten as a convolution of the injection spectrum dN^^/dEg in the 
channel with final state / with a halo function lie for the IG radiative process. We explicitly write 


d$!yC 

dDdE.y 


1 ^0 1 
2 47r ^7 



dEg Iic{E^, Eg, b, £) E( av) f 

f 


dN^ 

dEg 


(Eg) 


(4) 


where Eg is the electron (or positron) injection energy. In order to derive the differential y-ray flux produced 
by the secondary emission, we will compute Jic and dN^±/dEg, by using the tools in Ref. [l^ and then we will 
numerically perform the integral on the right-handed side of Eq. Q. 


^ The sum runs over the three different species of the ambient photon background; i.e. starlight, IR and CMB. For starlight and IR light, 
we extract the maps of their distribution from Galprop HZ). 
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FIG. 1: Several example of 7 -rays fluxes for 2 TeV DM candidates annihilating with thermal cross section ((ctv) = 3 x 
10“^®cm^/s) into leptonic channels (first row) and hadronic ones (second row). In all plots, the different hashing shows the 
predicted signal in two benchmark Rols with apertures 9i (dotted line), and O 5 (solid line). The red, blue and thick black lines 
denote instead the spectral features of the prompt, ICSs and total fluxes respectively. See the text for further details. 


B. Definition of the Rols relevant for our analysis 

Having at our disposal the fluxes per steradian, one needs to integrate Eqs. ([^ over a region AQ which cover, 
for instance, a given CTA observational window. In particular, we will consider five Rols corresponding to annuli 
centered in the Galactic Center with I-esime aperture Oi and constant thickness A0, minus rectangular regions with 
latitude within = 0° < b < &niax = 0-3° and i-esime longitude which intercepts the annulus thickness"^. The solid 
angle covered by the i-esime Rol will then be 

{ AH“” = 2n f d9 sin 9 (I-esime annulus) 

- (5) 

AH = 4 / d£ db cos b (i-esime rectangular region) 

and the corresponding averaged J-factor and averaged IC halo function Eg) = 

{Ej, Eg) — {Ej, Eg) are obtained by integrating the function J{9{b,£)) and Iic{Ey, Eg,b,£) over the 

region AH^ defined in Eq. ([^. For our morphological analysis, the annulus thickness is always A9 = 1° and we choose 


In the central 300 pc of the GC {i ~ 0°) and at Galactic latitudes with | 6 | < 0.3°, the H.E.S.S. collaboration discovered two bright 
sources: i) HESS J1745-290 coincident in position with the supermassive black hole Sgr A* |49| : and ii) HESS J1747-281 coincident 
with the supernova/pulsar wind nebula GO.9+0.1 |5Q| . together with a strong astrophysical diffuse emission |51| (see Sec. IV.B). This 
region is excluded in order to avoid VHE 7 -ray background contamination in the considered Rols. 
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five values of Oi rangine; from 1° to 5° with step of 1°. In the rest of the paper the label i = \,.... 5 will always refer 
to Rols with Oi = 1°,...,5°. 

Table |T] shows the angular size of the z-esime Rol together with the corresponding value of the averaged J-factor in 
units of GeV^/cm^. As one can see, although the angular size of the outer Rols is bigger, the values of the J-factor 
stay almost constant for regions above the second Rol. Averaged J-factors, assuming an Einasto profile with the 
same parameters adopted in Ref. [39] and a NEW profile as in Ref. |45|, are showed in the same table for a sake of 
comparison. As one can see, the Einasto profile used in this study provides conservative estimate of the J-factor in 
the considered Rols with respect to the same profile but with different parameters adopted in Ref. [31] . On the other 
hand, the NEW profile, which is cuspier at the GC than Einasto, provides a shightly smaller value of the J-factor in 
the considered Rol. This is mainly due to the fact that the Einasto profile is somewhat more chubby than NEW at 
few 100 pc to kpc scales (see e.g. Fig. 1 of Ref. [35]'). 

The integrated total flux (primary -f secondary) coming from the z-esime Rol will finally write 


d$tot ^ 

dE^ ’ 


In Fig. we show the y-ray fluxes for DM candidates with thermal cross-section ((CTu)th = 3 x 10“^® cm^/s) in 
leptonic channels (first row), and in hadronic ones (second row). In all plots the hashing shows the predicted signal in 
two benchmark Rols with apertures 9i (dotted line), and 9^ (solid line). Together with the values reported in Tab.[^ 
we can clearly see that the fluxes increase of roughly a factor 3 moving from the closest to the farthest Rols used in 
our morphological analysis. 

As an aside, the different plots in Fig. [l] allow us also to appreciate the different spectral features of the primary 
(red lines), secondary (blue lines) and total (thick black lines) emissions. As we can see, for a 2 TeV DM annihilating 
into e+e” and the y-rays flux spectra are dominated by the secondary emission up to a photon energy close 

to the DM mass. For the channel, the prompt emission is more pronounced and it dominates above a photon 

energy of roughly 100 GeV, while for the hadronic scenarios the secondary emissions are always subdominant in the 
entire CTA energy window. 

Finally it is worth noticing that for the channel DM DM— W~^W~, a pronounced spectral feature close to muM 
appears. This is due to the fact that the photons coming from the splitting —>■ IF^y contribute to the signal in a 
significant way. That electroweak process in fact, has two “soft” singularities: one comes from the usual soft photon, 
the other from a soft W (thus a y that carries away almost all the W energy). The latter divergence is screened by 
the mass of the W. Hence, we expect that the feature one can observe in the y-ray fluxes from the W^W~ channel, is 
only relevant for very heavy DM candidates (todm mw)- We stress here, that this pronounced spectral feature can 
be used for deriving constraints for lines-like searches has already done for the DM DM—yy or yZ by the H.E.S.S. 
collaboration jSj. 
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C. Expected photon counts in the CTA array 

In order to properly reproduce the expected photon counts in the CTA array we need to take into account the 
characteristics and performances of the detector. As stated in Sec. |ll|we consider the instrument response functions 
obtained for the proposed array I. For a given observation time Tobsj the number of observed y-ray events in the 
energy bins of the benchmark array is then obtained by convolving Eq. (§ with the energy-dependent effective area 
for photon and the Gaussian energy response. In the z-esime Rol and _)-esime energy bin with width AE^, 

it writes 




dE. 


lAEi 


dr® ■ 

7,2 

dE^ 
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7,2 

dE^ 


f-\-oo 
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'-{e')a:^{e') 
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where a-{E^) = S,-es{Ey)/ 78^(2) is the energy dependent Gaussian width and 7es(A.y) is the energy resolution taken 
from Ref. [13]. Here dr7/hE.y is the expected differential rate of counts. The left and central panels of Fig. 
shows the spectral morphology of E^ ■ d^^°\/dE^ and dF® ;^/dif.y respectively, considering 2 TeV DM candidates with 
thermal cross section into bb (black solid line) and e“'"e“ (red solid line) primary channels. From the comparison 
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First Rol, Einasto profile 


First Rol, Einasto profile 


Einasto profile 



ov = 3 X 10 cm s 
= 2 TeV 

— DMDM ^bb 

— DMDM -»e e 
Residual bck 
GDE 



E (TeV) 


0.1 0.2 1 2 3 45 10 20 
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• Residual bck X 10 

• GDEx lo" 


ov = 3 X 10 cm s 
m =2 TeV 


FIG. 2: Spectral and spatial morphology of typical DM signals together with the irreducible backgrounds. Left panel: Expected 
DM differential fluxes in the first Rol, considering 2 TeV thermal DM candidates annihilating into bb (black solid line) and 
e^e“ (red solid line) primary channels. For a sake of illustration, we show a spectrum for the CR hadrons (protons + nuclei) 
background (dashed green line) multiplied by a rejection factor e = 10“^ and a log-normal distribution for the electron spectrum 
(green dotted line) without rejection factor (ce = 1). The parametrizations of the irreducible backgrounds are given in Tab. 3 
of Ref. |42]. Notice that the energy E in the y-axis, refers to the “true” energy of the particle which initiates the shower. 
The GDE flux detected by Fermi-LAT extrapolated at high energy is also shown as dashed blue line. Central panel: Expected 
differential count rate in the first Rol considering 2 TeV DM candidates with thermal cross section into bb (black solid line) and 
e^e“ (red solid line) primary channels. The most up-to-date CR background from a CTA Monte Carlo simulation (solid green 
line) and the GDE (dashed blue line) detected by Fermi-LAT in the same region are also presented for the sake of comparison. 
Right panel: Expected total rate summed over all energy bins in the i-esime Rol and for the same benchmark DM candidates 
(red filled triangle for annihilating DM into e^e“ pairs, black filled square for bb primary channel). The irreducible CR hadrons 
background in the i-esime Rol is identified as green filled circles, while our “extreme” isotropic GDE as blue-filled stars. See 
the text for further details. 


between the left and central panels, we can see that, at low energy, the DM signals are suppressed due to the reduced 
CTA effective area. This effect is sizeable below roughly 100 GeV: indeed it can reduce the low-energy rate by about 2 
orders of magnitude. On the other hand, at high energy, the rates do not go to zero exactly at todm = 2 TeV. This is 
because of the finite energy resolution of the detector that subtracts part of the signal close to ttidm and redistributes 
it towards higher energy. 

The right panel of Fig. shows instead the spatial morphology of the total rates summed over all the energy bins 
above 30 GeV (T^ ^ = l/TobsE,^®„) in the i-esime Rol and for the same benchmark DM candidates (green filled 
triangle for annihilating DM into e+e“ pairs, red filled square for hh primary channel). It is worth noticing a sign 
change of the slope of the DM signal shape between the Rols 3 and 4 (see also Tab. |T]). We anticipate that the gentle 
spatial dependence of the DM signal with respect to the steeper one of the irreducible backgrounds (right panel of 
Fig.§ will allow for a substantial discrimination power that can be used to improve the CTA sensitivity. 


IV. CTA IRREDUCIBLE BACKGROUNDS 

In case of lACTs, the dominant backgrounds for DM searches at the Galactic Genter are the residual flux of GR 
hadrons and electrons and the GDE. In what follows, we discuss and quantify these irreducible backgrounds that 
erroneously pass the photon hardware trigger and analysis cuts. 


A. CR protons and electrons background events 

The numerous interactions between CR hadrons (protons -F nuclei) and the Earth’s atmosphere trigger hadronic 
showers which induce electromagnetic sub-cascades due to the decay of neutral pions into photons. The resulting 
electromagnetic cascade has a very different spatial morphology with respect to the one originated by CR electrons 
and photons, making possible a shower shape discrimination. Nevertheless, since the flux of CR hadrons is much 
larger than the one of CR electrons and photons, an electromagnetic cascade from a pion decay can be erroneously 
reconstructed as a 7 -ray if the rejection factor of the instrument is not very small. For a sake of illustration, the dashed 








green line in the left panel of Fig. shows the spectrum of the isotropic CR hadrons background taken from Ref. [H] 
in the first Rol, assuming a rejection factor e = 10“^ as done in Ref. [39]. As will become clear later on, assuming a 
constant e is not justified since it is energy dependent. Furthermore, it is also important to point out that the hadrons 
flux in the left panel of Fig. is shown as a function of the “true” energy of the particle which initiate the cascade. 
This is not the measured photon energy by CTA because of the reduced Cherenkov light emitted by the hadronic 
shower. Hence, an estimation of the CR hadrons background in the CTA energy window, can be obtained by shifting 
down the energy of the hadron. For example, assuming that only protons contribute to the CR hadrons flux, one can 
shift Ep to lower energies by a factor of 3, as done in Ref. [3S] (i.e. one can multiply the differential hadrons flux in the 
left panel of Fig. by a factor of 3 • 3“^'^ ~ 0.155), because among all the pions produced by the proton interactions 
with the Earth’s atmosphere, only 1/3 are neutral and in turn initiate electromagnetic sub-showers. Nevertheless, 
this rescaling is a rough approximation because also heavier species, especially He, are relevant (see e.g. Tab. 3 of 
Ref. |12|). Therefore, we stress here that the left panel of Fig. [2 can be only used for rough comparison. 

CR electrons constitute an annoying background for 7 -ray observations at Earth’s surface, because they initiate 
quasi indistinguishable electromagnetic cascades with respect to those originated by 7 -ray interactions. However, 
some discrimination between electrons and 7 -rays is possible by reconstruction of the primary interaction depth of 
the event jUj. For a sake of illustration, the dotted green line in the left panel of Fig. shows the spectrum of the 
isotropic CR electrons background taken from Ref. |42| in the first Rol, without rejection factor (ce = 1). 

The total number of residual background events (hadrons + electrons) is then given by: 
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where d<i>^^/(dE-ydH) is the total CR background flux per steradian and Ag^{Ej) is the energy-dependent effective 
area for CR. In our analysis, the counting rate of the irreducible background events dr^^/dE.^ is directly extracted 
from Ref. [42] and provides the most up-to-date background computation for CTA. In so doing, we neither assume 
the rejection factors for CR hadrons and electrons, nor the shift at low energy to take into account for the reduced 
Cherenkov light emitted by hadronic showers, because they directly come from a full CTA Monte Carlo simulation. 
The solid green line in the central panel of Fig. shows the spectral morphology of this rate in the first Rol for the 
considered array configuration. 

We have checked that our irreducible background, which comes from full CTA Monte Carlo simulations, is different 
compared to the analytic one implemented in Ref. [39] . To quantify the discrepancies, we can multiply the counting 
rate of the irreducible background events in the central panel of Fig. ^by the factor E^/{A^^{E^)Ani) and directly 
compare it with the solid black line in Fig. 2 of [331 ■ We get a higher xotal background especially above I TeV where 
the difference can be bigger than one order of magnitude. 

This improved background treatment will then degrade the sensitivity to DM annihilations compared to Ref. [39] . 
We think that such differences are due to three main reasons. First, the latter paper tried to model the irreducible 
background by introducing a constant rejection factor by hand on the incoming proton flux (in particular in Ref. [39] . 
ep = 10“^). This is a rough approximation because generally Cp depends on the energy. Second, they used the 
same effective areas for CRs and photons. This is an unjustified assumption, because the effective area depends on 
the particle which initiate the shower. Third, they assume that only protons are relevant. This is again not true 
because also the CR He interactions with the Earth’s atmosphere generate an important irreducible background in 
the CTA energy window. Hence, taking a constant efficiency, assuming that only protons are relevant and taking 
can lead to underestimate the total CR background. 

As an aside, the central panel of Fig. allows us also to appreciate the potential discrimination power from the 
peculiar spectral shape of the DM signal, characterised by either sharp energy cut-off or bump-like features, with 
respect to the smoother spectral shape of the residual background. The left and central panels of Fig. also shows 
the GDE 7 -rays (dashed blue line) background, which we discuss in detail in the next section. 

The right panel of the same figure shows instead the spatial morphology of the total rates of background events 
summed over all the energy bins above 30 GeV (F^^^ = i/TohsJ2j opposed to the DM signal, the 

background shape over the Rols (green filled circle) is monotonic. This is due to the fact that the background is 
roughly isotropic and therefore the expected background events increase proportionally with the size of the Rols 
= ADi/ADfc). From the different energy and spatial behaviours compared to the DM signals, we then 
foresee a substantial gain in sensitivity. 
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B. GDE background events 


The inner Galactic halo is a very crowded region with numerous astrophysical emitters in the VHE 7 -ray regime. In 
2006, besides the central 7 -ray emitter (HESS J1745-290 [13[52]) and the supernova/pulsar wind nebula GO.9-1-0.1 [50] . 
the H.E.S.S. experiment discovered a diffuse 7 -ray emission at energies of (0.2 — 20) TeV. This emission spreads through 
the central 200 pc along the Galactic plane, with a spatial extension of about ±0.3° in Galactic latitude, likely to 
be correlated with the Gentral Molecular Zone [5T]. In order to avoid contamination from the strong astrophysical 
background from the above-mentioned VHE 7 -ray emissions. Galactic latitudes with |5| < 0.3° are excluded in the 
Rols considered in this study as briefly stated in Sec. [HI B| 

On the other hand, at lower energies (below roughly 500 GeV), the Fermi-LAT satellite has measured a diffuse 7 -ray 
emission [531154] . originated mostly by the tt® decays from proton-proton collisions in the interstellar medium. In our 
Rols, the GDE is well described by a power-law spectrum extending up to roughly 500 GeV and an accurate mapping 
of this emission in the inner 5° can be in principle inferred from Ref. |55j . Nevertheless, since above the CTA threshold 
the 7 -ray contamination is not fully understood, we will limit ourselves by considering two extreme choices of the 
GDE background. In one case, we do not consider at all the GDE in the CTA energy window. In another case we 
will assume an isotropic flux all over the Rols using an averaged normalization for the GDE spectrum. In particular, 
we take the GDE flux per steradian d4>^°®/(dE.ydH) coming from the public P7V6 GDE model within one degree 
from the GG by the Fermi-LAT collaboration and we extrapolate it to higher energies to cover all the accessible energy 
range of CTA. This latter choice is conservative because we are taking the higher 7 -ray contamination in all Rols and 
therefore the CTA sensitivity to DM signals will be maximally deteriorated. Indeed, it is reasonable to expect that: 
(i) the GDE we are using is bigger than the expected one in the first Rol (its normalization is in fact taken from a 
region with the same angular aperture of the first Rol, without however eliminating the rectangular part which cuts 
the Galactic plane); (ii) the GDE decreases in the outer Rols (see e.g. Ref. |39j where, in their optimistic scenario, 
an exact knowledge is assumed of the spatial morphology of the GDE); (Hi) the GDE is completely absorbed at high 
energy due to photon-photon scattering and pair production on ambient photon background radiation (e.g. the Fermi- 
LAT collaboration has recently reported a possible energy cut-off around 800 GeV, in the spectrum of the isotropic 
diffuse emission [5T]T 

Having at our disposal the flux, the number of GDE background events is again given by 
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where dT^^^/dE^ is the differential count rate. The blue dotted lines in the left and central panels of Fig. show 
the spectral morphology of E^ ■ d$®d®"/dE.y and dr^d®/dE.y, respectively. From the central panel, we clearly see 
that the GDE will increases the total background up to a few TeV energies. This will impact our expected sensitivity 
especially for low DM masses as we will show in Sec. |VI| 

The filled-blue stars in the right panel of Fig.[^ shows instead the spatial morphology of the total GDE rate summed 
over all energy bins above 30 GeV (F®^® = ^/TohsJ^j Since, we assume that the GDE is isotropic, the 

spatial dependence of this background is clearly contained in AD^. 


V. ANALYSIS METHODOLOGY 

As stated in Sec. [Hj observations by lAGTs in high energy 7 -ray astronomy usually employ two regions on the sky: 
one in the direction where the signal is expected to be (ON region), and one for the background determination (OFF 
region). As follows, we will implement the Li&Ma notations [56] used for the above mentioned ON-OFF method, by 
carrying out a full likelihood analysis which uses the expected spectral and spatial morphology of the DM signals. 


A. Sensitivity calculation methodology 

The statistical method to derive the sensitivity is based on a likelihood ratio statistical test. In order to take full 
advantage of the specihc spatial and spectral features in the DM signal (see Figs. ( ]I|2[ )), the total likelihood for a given 
todm is obtained from the product over the spatial bins i and the energy bins j of the individual Poisson likelihoods. 
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FIG. 3: CTA sensitivity to DM annihilation in the (moM, (rrn)) plane. Unless otherwise stated, our morphological analysis 
employs five adjacent Rols, the exposure in each Rols is 500 h, the energy threshold is 30 GeV and only the statistical 
uncertainties are taken into account. Upper-left panel: Improvement of the CTA sensitivity for the DM DM —>■ bb channel due 
to our multi-bin morphological method (two, three, five and eight Rols) compared to the case where only the energy spectral 
information either in the first (solid black line) or in the second (dashed black line) Rol is taken into account. Upper-right 
panel: CTA sensitivity for the DM DM —> bb channel when only the spectral information in the individual Rols (Rol 1 to Rol 
5) is considered. For a sake of comparison, we also show how our multi-bin morphological method (solid red line) ameliorates 
the CTA sensitivity. Bottom-left panel: Degradation of the CTA sensitivity for annihilating DM into bb pairs due to possible 
systematic errors in the rich observational datasets. Observation times of 100 h (black lines) and 500 h (red lines) assuming 
different values of systematics: 0.3% (dashed lines), 3% (dotted lines), only statistical fluctuations (solid lines) are taken into 
account. Bottom-right panel: Impact of an “extreme” GDE added on top of the GR background, on the CTA sensitivity. The 
sensitivity is shown for annihilating DM into bb (back lines) and e^e“ (red lines) pairs with (dotted lines) and without (solid 
lines) the irreducible GDE background. See the text for further details. 
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It explicitly writes 


CirriT^M, {o-v}) = ( ctu )) , 


( 10 ) 


* j 


where the spectral part of the likelihood runs from the energy threshold of 30 GeV up to ruDM, while the spatial part 
runs over the five Rols described in Sec. |IIIB[ Following Ref. [57] . the individual likelihood is then given by 


(-^7 ) -^7 l-^ON, -^OFf) = 




Nqff,'. 




( 11 ) 


where is the predicted number of events in the i-esime Rol and j-esime energy bin which corresponds 

to the observed counts in the ON region IVoN.y ■ Here is the predicted number of background 

events which is estimated from an OFF region in order to avoid modelling of the background contamination in 
the considered Rol. The observed number of background events in the OFF region is denoted by NoFF,ij and the 
parameter ai = AOi/AOoFF refers to the ratio between the angular size of the Fesime ON region and the OFF one. 


From Eq. (111, one can then infer that NoN,ij and fVoFF.ij correspond to Poisson realizations with mean ^ ■ + A® 


and N^^/ai respectively. The background is estimated beyond the fifth Rol in a region five times larger. In this 
case the ratio between the fifth Rol and OFF region is Ui = 0.2. We assume a conservative constant value of a even 
if the value of it is smaller for the inner four Rols. For our statistical analysis we will adopt the likelihood ratio test 
statistic TS = —2 ln(£(TODM) (o'u))/£niax(™DM, (o'v))) which follows an approximate distribution with one degree 
of freedom. Values of TS higher than 2.71 are excluded at a 95% Confidence Level (C.L.). 

Another possible approach, which has been used e.g. in Ref. [31], is to consider only ON regions in the likelihood. 
In this case, the Poissonian term for Aoff.ij in the right-hand side of Eq. ( [II| ), is not considered at all. Therefore 
the background determination relies on a careful modelling in the Rols. Although this approach allows to put 
stronger constraints, it crucially depends on the residual background modelling accuracy. In our multi-RoI approach, 
the background is computed from data taken in real observation conditions and therefore allows for more accurate 
background determination. 


B. Impact of background systematics 


CTA observations towards the inner Galactic halo will provide statistically rich observational datasets. This datasets 
result from the combination of multiple observations with distinct observational parameters that may introduce specific 
observational systematics. 

Since the systematic errors may be the limiting factor in the overall error budget for GTA, an assessment of the 
impact of them in the GTA sensitivity is in order. A possible way to investigate the impact of such uncertainties is 
to introduce in the likelihood a Gaussian nuisance parameters m such as the individual likelihood writes 




, /3\Nof!, Noff) = 
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( 12 ) 

where /3ij acts as a normalization parameter and (Ty is the width of the Gaussian function. An accurate determination 
of the spatial and energy dependencies of the systematic level is beyond the scope of this study, and we will keep Uij 
fixed for all spatial and energy bins. We can determine the maximum likelihood value of it by solving AC/df^ij = 0 in 
the z-esime Rol and j-esime energy bin, for a given set of A® AoN.ij and Aoff.ij- In order to evaluate the 

impact of systematics in the GTA sensitivity, we will then consider several plausible values of the Gaussian width a 
for lAGT observational data. In particular we will vary it from 0.3% to 3% in order to compare with the case where 
only the statistical uncertainty is taken into account. 


VI. RESULTS 

We now show our results in terms of 95% G.L. sensitivity limits on DM annihilation, in the usual (ttidm, ctu) plane. 
We focus on several particle-antiparticle annihilation modes (DM DM —>■ e'^e~, ,bb,tt and W'^W~) and 

in a broad range of DM masses (from 30 GeV up to 80 TeV). For all the channels, we include the ICSs of energetic 
produced by annihilating DM on the ambient photon background which is particularly relevant for e“'"e“ and 
modes (see the red lines in Fig.[^. We summarize our main results in Fig. In particular we find that: 
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o The ICS emission substantially increases the CTA sensitivity for the leptonic channels. In Fig. we show the 
spectral features of the ICS emissions. As it is apparent, the total fluxes receive a substantial contribution for 
photon energy just below the DM mass in case of leptonic channels (especially for the DM DM —>■ e+e” and 
pCp,~ modes). As a consequence, since the IC secondary emission is well inside the CTA energy window, the 
sensitivity to those channels is largely ameliorated. 

o A morphological analysis ameliorates the CTA sensitivity. Assuming again an observation time of 500 h, in the 
upper-left panel of Fig. we show the improvement of the sensitivity for the bb channel, by using the multi-bin 
morphological method compared to the case where only the energy spectral information either in the first Rol 
(solid black line) or in the second one (dashed black line) is used. The constraints inferred with this approach 
are more stringent by a factor of a few with respect to those obtained considering only the spectral information 
in any single Rol. Indeed, from the upper-right panel of the same figure, one can quantify the impact of our 
morphological analysis by comparing the different black lines (one Rol, Rol 2, 3, 4, 5 alone) with the red one 
(combination of all the relevant Rols considered in our analysis). As one can see all the constraints coming from 
individual region, are a factor 0(2) less stringent than the one inferred with a multi-bin morphological method. 
As shown in the right panel of Fig. this is mainly due to the fact that the spatial dependence of the DM 
signal flatten above roughly the third Rol, while the CR background contamination increases with the size of 
the f-esime region. 

o The systematic uncertainties deteriorate the CTA sensitivity. In the bottom-left panel of Fig. we show the 
impact on the CTA sensitivity for annihilating DM into bb pairs due to possible systematic errors in the rich 
observational datasets. We use our morphological analysis (hve Rols) with observation times of 100 h (black 
lines) and 500 h (red lines) assuming different values of systematics: 0.3% (dashed lines), 1% (red dotted line), 
3% (black dotted line), only statistical fluctuations (solid lines). We can see that the sensitivity is deteriorated 
over all the mass range for a given observational time. In particular, introducing a systematic errors of 0.3%(3%) 
for 100 hours, deteriorates the sensitivity of a factor 1.5 (see Fig. [^. For multi-TeV DM masses, the impact is 
reduced because at higher energy the y-ray datasets are dominated by the statistical errors. We stress that the 
systematics must be controlled at the level of 0.3% or better to allow for a substantial sensitivity improvement 
with a 500 h observation time. 

o Our extreme choice of the Calactic Diffuse Emission degrades the CTA sensitivity. Assuming an observation 
time of 500 h, in the bottom-right panel of Fig. we assess the CTA sensitivity for annihilating DM into bb 
(back lines) and e“''e“ (red lines) pairs once an “extreme” GDE is added to the CR background. As commented 
upon in Sec. |IVB| we consider an isotropic GDE coming from the inner Galactic halo. We can see that the CTA 
sensitivity is significantly deteriorated below DM masses of few TeV (dashed lines) with respect to the scenario 
where the GDE is not considered at all (solid lines). On a more specific level, even with our “extreme” choice 
of the GDE, the CTA sensitivity still probe cross section below the thermal value for the e+e” channel. For 
the hadronic ones (e.g. bb mode), the CTA sensitivity is degraded of a factor 2 making the reach of the thermal 
cross section no longer possible. Nevertheless, since we assume that the GDE is isotropic, it is worth stressing 
once again that we are overestimating the y-ray contamination in all Rols. In fact, if we consider an accurate 
mapping of the GDE in our Rols (like the one used in Ref. |39j in their optimistic scenario), we find that the 
impact of the GDE in the final results is not very pronounced. This is due to the fact that in all the regions used 
in our analysis, the diffuse y-rays contamination is smaller than the residual CR background extracted from a 
full CTA Monte Carlo simulation. 

In Fig. unless otherwise stated, we show the most optimistic limits in the (ru-DMi (cu)) plane in which neither the 
systematic uncertainties in the datasets nor the GDE are taken into account. 

In particular, in the left panel of Fig. |^we show the CTA sensitivity to DM annihilation assuming different channels 
(DM DM —)■ e'^e~, pcfjT,T^T~ ,hb,ti and W'^W~), and an observation time of 500 h. Focussing first on the purely 
leptonic channels, we find that CTA would be able to exclude annihilation cross-sections well below the thermal value. 
On a more specific level, the best sensitivity is obtained for the DM DM —e+e“ mode {{av) ^ 5 x 10“^^ cm^/s 
for todm — 200 GeV) rather than pCpT and r+r“, since the original population is produced at higher energies, 
and therefore the secondary ICS emission is well inside the CTA energy window (see the upper raw in Fig. [^. For 
the hadronic and W~^W~ channels, we get the same qualitative feature of the exclusion limits modulo a factor of 
©(few) in the normalization. This can be explained from the fact the y-ray spectrum arising from the fragmentation 
of sufficiently heavy hadronic SM particles is quasi-universal. In particular, for the DM DM —>■ bb channel the 
thermal value of the cross section can be probed in the TeV mass range, where the best sensitivity is achieved at 
(av) ~ 2 X 10“^® cm^/s for ttidm — 1 TeV. 
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Five Rols, Einasto profile, 500 h 




FIG. 4: Unless otherwise stated, the most optimist limits in the (rriDM, {o'v)) plane in which neither the systematic uncertainties 
in the datasets nor the GDE are taken into account. Left panel: CTA sensitivity to DM annihilation for various primary channels 
(DM DM —>■ bb in black, it in blue, W^W~ in yellow, t'^t~ in green, pT in red and e''"e“ in dashed black) together with the 
reference value of the thermal cross section (long-dashed black line). The sensitivity is computed for a 500 h fiat exposure over 
the five Rols used in our analysis assuming a 100% branching ratio in each annihilation channel. Right panel: CTA sensitivity 
to DM annihilation for the DM DM —> bb channel (black line) compared to the most constraining limits to date. In particular, 
we report the H.E.S.S. limit for 112h of observations of the GC region (solid red line) and the Fermi-LAT limit coming from 
5 years of observations of 15 dSphs (solid blue line). An estimate of the projected sensitivity of Fermi-LAT for 45 dSphs and 
10 years of observation time is also shown (dotted blue line). For a sake of comparison, we also report the CTA sensitivity 
assuming a reasonable 0.3% of systematics (dashed black line). See the text for further details. 


We comment here on the relative strength of our constraints for the DM DM —>■ bb channel in Figs. ( |3|4[ ), with 
respect to the best limits to date obtained from either other analyses or targets. For a sake of comparison, in the right 
panel of Fig. we also report the CTA sensitivity assuming a reasonable 0.3% of systematics (dashed black line). 

With respect to the H.E.S.S. limits from 112 h of observations in the central region of the Galactic halo [7] (red line 
in the right panel of Fig. |^, the bounds derived here allows to probe a larger range of DM masses (from 50 GeV up 
to 80 TeV). In particular the CTA sensitivity with no GDE and systematics will be roughly a factor 10 more sensitive 
for DM masses around 1 TeV. 

With respect to the most recent stacking analysis on 15 dwarf Spheroidal galaxy (dSph) observations from Fermi- 
LAT (solid blue line in the right panel of Fig. |^, CTA becomes competitive for energies above 100 GeV and it 
overtakes the Fermi-LAT constraint for DM masses above 200 GeV. Furthermore, it is important to point out that 
CTA will still provide stronger limits above roughly 500 GeV, if the very optimistic scenario of 45 observed dSphs for 
10 years Fermi-LAT observations [SS] will be considered (dashed blue line in the right panel of Fig. [^. Hence, together 
with the optimistic Fermi-LAT observations of dSphs, CTA will be able to survey thermal DM candidates in a broad 
range of masses (from few tens of GeV up to several tens of TeV). 

With respect to previous CTA projected limits from Refs. |35[l37l[55] . our results are more conservative. In particular, 
in contrast with these works, we include both the systematic uncertainties in the datasets and an “extreme” GDE on 
top of the GR background. We find that the shape of the bounds as a function of the DM mass is very different. This 
is mainly due to the fact that we use both spectral and spatial information in the likelihood. The impact of including 
spectral information was pointed out in Ref. |38j . whose limits have in fact a similar shape compare to our bounds. 
The fact that our limits with no GDE and systematics are instead less stringent than those obtained in Ref. |35], is 
mainly due to the fact that we use a less steep profile rather than a profile inferred from the Aquarius simulation 
which include a significant boost due to substructure. Indeed, in their ON region, which almost coincide with our 
first Rol, the value of the J-factor based on the Aquarius simulation is 4.68 x 10^^ GeV^/cm®, while in our case is 
around 1.42 x 10^^ GeV^/cm®. This yields a difference in the flux compare to our work of a factor ~ 33. Nevertheless, 
since we use a morphological analysis which improves the CTA sensitivity of roughly a factor 2, we should multiply 
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the limits in Ref. [35] by 33/2 for a sake of comparison. 

With respect to Ref. |39j . which is similar to our work in many aspects, we find again that our results are more 
conservative. This is mainly due to the fact that the parameters of the Einasto profile used in our analysis are different. 
Considering their J-factor (in Ref. [39] the J-factor in our first Rol® (see Tab. is around 4.61 x 10^^ GeV^/cm®), 
and the fact that we implement, in a different way, a morphological analysis in the likelihood (improvement of the 
CTA sensitivity of a factor 2), we find however that the bound for the DM DM —bb channel for 100 h observational 
time is in well agreement with our own. The only notable difference has been found for DM masses above roughly 1 
TeV, since, as pointed out many times, we use the most up-to-date CR background obtained from a full CTA Monte 
Carlo simulation. This gives a larger contamination above roughly 1 TeV compare to the analytic estimation of the 
CR background used in Ref. [33]. As a consequence, our bounds are less stringent. The fact that our background is 
bigger at larger energies is also relevant for the discussion related to the GDE. In fact we find that only if we choose 
an “extreme” GDE, the constraints are significantly affected. If we use the accurate mapping of the GDE used in 
Ref. |39j in their optimistic scenario^ we only report a barely modification of the CTA sensitivity. Finally, concerning 
the impact of the systematics errors on the constraints, we find an almost perfect agreement with that presented in 
Ref. lag. 


VII. SUMMARY 

We discuss the future CTA sensitivity to DM annihilations in several channels and over a range of DM masses 
from 50 GeV up to 80 TeV. Eor all channels we include the ICS emissions which yield a substantial contribution 
to the overall 7 -ray flux, especially for the leptonic channels. We improved the analysis over previous work by: i) 
implementing a spectral and morphological analysis in the 7 -rays emission; ii) taking into account the most up-to-date 
CR background obtained from a full CTA Monte Carlo simulation and a description of the GDE on top of this; and 
Hi) including the systematic errors in the rich observational CTA datasets. 

We showed that our morphological analysis with five Rols substantially improves the CTA sensitivity by roughly 
a factor of 2. In particular, for the hadronic channels, we found that the CTA with an uniform exposure of 500 h of 
observations, will be able to probe thermal values of the annihilation cross-section over a broad range of DM masses, if 
the systematics uncertainties in the datasets will be controlled better than the percent level. For the leptonic channels, 
and in particular for the DM DM —>■ e+e” mode, we found that the bounds are instead well below the thermal value 
of the annihilation cross section. In this case, even with larger systematics, thermal DM candidates up to masses of 
few TeV will be easily probed. 

In contrast with Ref. [3g, we found that the inclusion of the GDE in the overall background produces a sizeable 
effect on the CTA sensitivity, only if one assumes an isotropic 7 -rays contamination extracted from the inner galactic 
halo. With a more physical choice, which includes a proper spatial dependence of the GDE, we only found a minimal 
impact in the CTA sensitivity. This is due to the fact that our CR background is larger than the one analytically 
inferred in Ref. [31] . Therefore, apart from an “extreme” choice of the GDE, the degradation of the CTA sensitivity, 
especially at high energy (above roughly 1 TeV), is mostly due to the CR background contamination. 

In summary, in order to probe the thermal value of the annihilation cross-section over a broad range of DM masses, 
deep observations of the GC over several degrees in radius (at least up to 5 degrees) are required with uniform exposure 
and the best possible control of systematic uncertainties. Under these conditions, CTA will give crucial information 
for TeV-ish WIMP searches in the next decade. Furthermore, in the optimistic scenario where the LHC will discover 
new physics in which thermal DM candidates are present, CTA will be probably the only player that could cross-check 
such possible results against an astrophysical environment. 
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^ Notice that in Refs. m the value of the J-factor is quoted in a region slightly bigger than ours. In particular, they considered a circle 
centered in the GC with an aperture of 1.3°, minus a rectangular region which cuts the galactic plane. 
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